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Abstract 

The effect of additive white noise on a model for bursting behavior in large aspect- 
ratio binary fluid convection is considered. Such bursts are present in systems with 
nearly square symmetry and are the result of heteroclinic cycles involving infinite 
amplitude states created when the square symmetry is broken. A combination of 
numerical results and analytical arguments show how even a very small amount 
of noise can have a very large effect on the amplitudes of successive bursts. Large 
enough noise can also affect the physical manifestations of the bursts. Finally, it is 
shown that related bursts may occur when white noise is added to the normal form 
equations for the Hopf bifurcation with exact square symmetry. 
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1 Introduction 



In physical systems, noise is unavoidable. Thus, it is important to study the 
effect of noise on models of physical systems. It has been shown that some 
models are very sensitive to the presence of noise: for example, the model for 
the resonant interaction of three wave modes considered in [1] shows that a 
tiny amount of noise can replace a bifurcation structure involving a period 
doubling cascade to chaos with a noisy periodic orbit which is attracting over 
a wide range of parameter values. Noise can have an especially profound effect 
on models with structurally stable heteroclinic cycles. These cycles involve 
connections between unstable states, and the connections typically occur in 
invariant subspaces. Even small amounts of noise can cause the system to 
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jump across these subspaces, and this can lead to very different dynamical 
behavior [2] . If such a heteroclinic cycle also is attracting, the time spent in the 
neighborhood of the unstable states of the cycle increases without bound with 
increasing time. This is not observed in physical systems and noise provides 
an explanation: the heteroclinic cycle becomes a "statistical limit cycle" in 
which switching between the unstable states of the cycle occurs randomly in 
time but with a well-defined mean period [3-5]. 

In [6,7], a model of large aspect-ratio binary fluid convection was studied. 
The model considers the competition between two nearly degenerate modes 
of opposite parity. The evolution equations correspond to a Hopf bifurca- 
tion with broken D 4 symmetry, where the "interchange" symmetry between 
the two modes is weakly broken because of the large but finite aspect-ratio. 
For open parameter regimes, periodic or irregular bursting with very large 
dynamic range can occur close to threshold. The bursts are associated with 
nonattracting heteroclinic cycles involving solutions "at infinity" [7,8]. These 
cycles may be structurally stable or unstable. Bursts occur when the sys- 
tem evolves toward an infinite amplitude state along its stable manifold, then 
gets kicked toward another infinite amplitude state with an unstable mani- 
fold which returns the system to finite amplitude. This mechanism provides 
an explanation for bursts found experimentally for large aspect-ratio binary 
fluid convection [9], and may also provide a deterministic explanation for the 
variability of the solar magnetic cycle [10]. A comparison of this to other mech- 
anisms for bursting is given in [11,12]. Related bursts can occur due to the 
resonant temporal forcing of a system undergoing a Hopf bifurcation with D 4 
symmetry [13]. 

When noise is added to the model considered in [6,7], there will be a burst with 
larger amplitude if the noise kicks the system closer to the stable manifold of an 
infinite amplitude state, and a burst with smaller amplitude if the noise kicks 
it away. Thus a state which bursts in the absence of noise is still expected 
to burst in the presence of noise, but the successive burst amplitudes and 
the time between successive bursts may change substantially. The details of 
the effect of noise on the bursting behavior are described in this paper. This 
includes the possibility that noise can affect the physical manifestations of 
the bursts. Specifically, purely blinking and winking states in which successive 
bursts occur at opposite sides or the same side of the container, respectively, 
can be destroyed by sufficiently large noise. Finally, it is shown that related 
bursts can occur even if the D4 symmetry of the model is not broken, i.e., 
when noise is added to the normal form equations for the Hopf bifurcation 
with exact D 4 symmetry. 
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2 Evolution equations for the model 



For binary fluid convection, if the separation ratio of the mixture is sufficiently 
negative then convection arises via a Hopf bifurcation. This is the case for the 
3 He/ 4 He mixture used by Sullivan and Ahlers [9] in their experiment carried 
out in the large aspect-ratio rectangular container 

D = Ix, Y, Z\ — \v < X < \t, -\t y <Y < \v Y , -\ < Z < \\ , 
1 2 ~ - 2 ' 2 - 2 ' 2 _ - 2J ' 



with aspect-ratio T = 34, and Ty = 6.9. In this experiment, it was observed 
that immediately above threshold convective heat transport may take place 
in a sequence of irregular bursts of large dynamic range despite constant heat 
input. A model for this experiment is considered in [6,7]. The perturbation to 
the conduction state temperature profile is assumed to take the form 

G(X, Y, Z, t) = e^ 2 Re{z + (t)f + (X, Y, Z) + z_(t)f_(X, Y, Z)} + (9(e), 



where eC 1, f±(—X, Y, Z) = ±f±(X, Y, Z). Here z + are z_ are the (complex) 
amplitudes of the first modes to lose stability which are respectively even 
and odd under the reflection X — > —X. In [6], evolution equations for z + 
and z_ were derived using symmetry arguments, with the resulting equations 
describing a Hopf bifurcation with broken D4 symmetry. The applicability of 
this model to the experiment is discussed further in [14]. As shown in [6,7], 
solutions of the evolution equations include bursts of very large dynamic range 
similar to those found in the experiments. The mechanism by which these 
bursts arise is discussed in detail in [8], and is summarized below. 

In this paper we generalize the model considered in [6,7] to include additive 
random forcing terms. Specifically, we consider the formal equations 



(A + AA + i(u + Au))z + + A(\z + \ 2 + \z^\ 2 )z + 
+B\z + \ 2 z + + Cz + z 2 _+ri 1 (t)+iri 2 (t) (1) 
(A - AA + i(u - Au))z_ + A(\z + \ 2 + \z-\ 2 )z^ 
+B\z_\ 2 z- +Cz^z 2 + +r] 3 (t) +ir] 4 (t) (2) 

where the r/^s represent real, independent Gaussian white noise random pro- 
cesses with the properties 

( Vl (t)) = 0, ( Vi (t)r 1j (t')) = 2V8(t-t')5 ij , i,j = 1,2,3,4. (3) 



dz + 
~~dt 



dz_ 
~dT 
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Here z+, Z-, A, B, C are complex, and A, AA, u, Au, T> are real. The quantity 
Au measures the difference in frequency between the two modes at onset, and 
AA measures the difference in their linear growth rates. The terms involving 
AA and Au are called forced symmetry-breaking terms because they break 
the D 4 symmetry of the governing equations [7]. The new rji terms represent 
unavoidable random effects in the experiment. 

We first summarize results for the case that there is no noise, i.e., when T> = 0. 
New coordinates (p, 9, 0, ip) are defined according to 

z+ = p -V2 cos (^/2) e ^)/ 2 , = p- 1 ' 2 sin(£/2)e l (-^/ 2 , 

where without loss of generality 9 G [0,7r], G [— 2tt, 2tc), and -0 G [0, An). 
Also introducing a new time r defined by dr/dt = 1/p, equations (1,2) with 
T> = take the form 

-p[2A R + B R (1 + cos 2 9) + C R sin 2 9 cos 2(f)] - 2(A + AA cos 9)p 2 (A) 
sin #[cos 9{-B R + C R cos 20) - d sin 20] - 2AA sin Op, (5) 
cos 9(Bj - Cj cos 20) - C R sin 20 + 2Acup, (6) 
2Ai + Bi + Ci cos 20 + C R sin 20 cos 9 + 2up. (7) 
A# + etc. Let 

1 I 1 2 i I 1 2 

r = - = z + + z_ 
P 

denote the amplitude of a solution. As discussed in [7,8], equations (4-7) have 
an invariant subspace X at p = (corresponding to infinite amplitude states) 
on which the equations are equivalent to those with A A = Au = 0. Since 
the variable ip decouples from the other variables, the dynamics on S are 
two-dimensional and hence simple to analyze [15]. Such an analysis (see [8]) 
allows us to conclude that, for this truncation, infinite amplitude periodic and 
quasiperiodic solutions exist. Bursts are associated with heteroclinic cycles in- 
volving such infinite amplitude states. Reference [8] shows that for such cycles 
to form, it is necessary to have at least one subcritical and one supercritical 
solution branch bifurcating from the trivial state with AA = Au = 0; also, the 
"angular stability" properties of the solutions must be correct. Reference [8] 
also shows that despite their heteroclinic nature, the duration of the resulting 
bursts in the original time t is in fact finite, i.e., the excursion to and return 
from infinity occur in finite time. Of course, the infinite amplitude solutions 



dp 

d^ = 
d9__ 

d^~ 

dcf)_ 

dr 
dr 

Here A = 
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and connections to them are of physical interest only insofar as they are re- 
sponsible for the presence of nearby solutions which make visits to large but 
finite amplitude. This type of bursting behavior persists for T> = even when 
higher order terms in equations (1,2) are retained [8]. It should be emphasized 
that the noise is added on the original timescale t, not r (see equations (1,2)). 



3 Effect of noise on bursting behavior 

We now describe the effect of noise on the bursting behavior reported in [6,7]. 
Coefficient values which give the appropriate sub- and supercriticality of trav- 
eling and standing wave branches in large aspect-ratio binary fluid convection 
systems are considered [7]: 

A = l-1.5i, B = -2.8 + 5i, C=l+i, 

co = l, AA = 0.03, Au = 0.02. 

For these coefficient values in the absence of noise, nonattracting, structurally 
stable heteroclinic cycles connecting infinite amplitude periodic orbits (here- 
after, Uoo solutions) and infinite amplitude quasiperiodic orbits (hereafter, 
qpoo solutions) exist over a range of A values. These heteroclinic cycles are of 
Shil'nikov-Hopf type, and there are associated stable solutions which display 
bursting behavior [7,8]. 

For definiteness, consider the stable quasiperiodic solution which is present 
in the absence of noise for A = 0.1 (see Figure 1(a); in the notation of [8], 
this is a w\ solution). The trajectory makes successive visits near different 
but symmetry-related solutions given by (p, 9, 0, ip) = (0, ir/2, inn, (2Ai + 
Bj + Ci)t mod Atc), m = 0,±1,±2. The Uoo solutions are unstable within 
S, and the trajectory gets kicked toward a qp^, solution. The qp^ solutions 
are unstable in the p direction, so the trajectory returns to smaller amplitude 
(larger p). Physically, this solution corresponds to a blinking state for the 
convection system because successive bursts occur at opposite sides of the 
container [7]. Figure 1(b) shows the corresponding results for V — 1 x 10~ 7 . 
Clearly, even this very small amount of noise can have a very large effect 
on the burst amplitudes. These and other numerical results were obtained 
using a stochastic second-order Runge-Kutta algorithm [16]. Unless otherwise 
indicated, the time step of integration is 5t = 1 x 10~ 4 . 

The results shown in Figure 1(b) may be qualitatively understood by recog- 
nizing that there will be a burst with larger amplitude if the noise kicks the 
trajectory closer to the stable manifold of a solution, and a burst with 
smaller amplitude if the noise kicks it away. This is elucidated by ignoring the 
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Fig. 1. Bursts for A = 0.1 and (a) T> = 0, (b) T> = 1 x 1CP 7 . Clearly, even a very small 
amount of noise can have a very large effect on the burst amplitudes. Fixed points 
(periodic orbits) in these projections correspond to periodic orbits (quasiper iodic 
orbits) for the full system. In this and later figures, the Uoo solutions are denoted 
by +'s, and the qpoo solutions appear as periodic orbits at p = 0. 

uncoupled variable tp, linearizing equations (4-6) about the solutions, and 
denning coordinates (p, x, y) to give the diagonal form 

dx dp dy 

dr = dr = P ' dr = ~ ^ 



where \ u ,\ s ,\ ss > 0. The stable manifold of the solutions then ap- 
proximately intersects a plane of sufficiently small, constant p along the line 
x = 0. For the parameters under consideration, A M = 0.0682, \ s = 0.2, and 
\ ss = 5.868, and equation (8) is obtained by using the coordinates approxi- 
mately given by 



x 



1 

0.269 -0.730 0.706 
0.00626 0.902 0.467 



P 

o-ti/2 

6 — 17171 



\ 



Figure 2(a) shows values at which the trajectory intersects the surface p — 1 
with decreasing p for T> — 1 x 10~ 7 . This is interpreted as showing how the 
noise affects the proximity of the trajectory to the stable manifold of the 
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Fig. 2. (a) Values of (x, y) at which the trajectory intersects the surface p = 1 with 
decreasing p for T> = 1 x 10~ 7 . The stable manifold of the Uqq solutions intersects 
the surface p = 1 approximately along the line x = 0. (b) Probability distribution 
function p(x) at the surface p = 1. The histogram is from a long-time integration of 
equations (1,2), and the solid line shows a fit to equation (9). (c) Variance a 2 from 
fitting p(x) to equation (9) for different noise strengths V. 

solutions. The probability distribution function (pdf ) of x values at the surface 
p = 1 is shown in Figure 2(b). By definition, p(x)dx is the probability that 
the x value at the surface p = 1 is between x and x + dx. This is fit reasonably 
well by the Gaussian distribution 

p(x) = . 1 exp \-(x - x) 2 /(2al)] . (9) 



Figure 2(c) shows that the variances a 2 x of fitted Gaussians vary linearly with 
the noise strength V. Over this range x 0.201, independent of V. The noise 
strengths have been restricted to the range shown in the figure because for 
larger T> the noise can kick the trajectory close enough to the stable manifold 
of a solution that the long-time integration of equations (1,2) becomes 
numerically difficult. 

We now show how to understand the statistical properties of the bursts in 
terms of the properties of p(x). In the following, we ignore the strongly con- 
tracting y direction. Define S p = {(x, p)\p = p*} and S x = {(x,p)\x = x*}, 
and consider how the pdf p(x) at S p evolves into a pdf p(p) at T, x (see Fig- 
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4- 



Fig. 3. Setup for determining how the pdf of x values at constant p = p* evolves 
into a pdf of p values at constant x = x* . 

ure 3). When the amplitude r becomes large (i.e. when p becomes small), 
the noise terms in equations (1,2) are overwhelmed by the cubic deterministic 
terms. Therefore, provided p* is sufficiently small, we can ignore the noise over 
the part of the trajectory shown in Figure 3 and approximate the vector field 
by equation (8). The point (xq, p*) G S p is mapped into (x*,f(xo)) G Y, x un- 
der the flow, where f(x) = p*(x/x*) Xa ^ Xu . Given p(x), the pdf p(p) is obtained 
using \p(x)dx\ = \p(p)dp\, with p = f(x). For sufficiently small noise we can 
take f(x) Ki f(x) + f'(x)(x — x), where x is the average value of x at T, p . 
Using (9) as the pdf for x at E p , the pdf for p at Y, x is then 



P{p) 



: CXp 



-(p-p) 2 /(2<7 2 ) 



In particular, p(p) is Gaussian. 



(10) 
(11) 



Now it is assumed that in the trajectory's successive visits near solutions, p 
always reaches its minimum value p min at some fixed x = x*. This is reasonable 
because at some x value the trajectory will have been kicked close to a qp^ 
solution which is unstable in the p direction, so the trajectory starts to evolve 
towards larger p (cf. Figure 1). With this assumption and from (10), the pdf 
for p min is 



p(Pmin) = 



\phxo'- 



: cxp 



-(Pmin - Pmin) 2 /(2CT 2 ) 



(12) 



The pdf for the maximum burst amplitudes, r max , is then obtained using 

^max = 1/Pmin and |p(r max )rfr max | = |p(p min )dp min |: 
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P(^max) 




Fig. 4. Probability distribution functions for (a) r max , (b) p m m, and (c) T for 
A = 0.1, P = 1 x 10~ 7 . The histograms are from a long-time integration of equa- 
tions (1,2), and the solid lines show fits given by equations (12-14). 
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Fig. 5. Linear dependence of the variances (a) a 2 , (b) a\ on the noise strength V. 



This distribution has a long 1/V4 ax tail for large r max , indicating the probabil- 
ity of bursts having very large amplitudes. These functional forms for the pdf's 
are verified for V — 1 x 1CT 7 in Figure 4(a,b). For the range of noise strengths 
considered, p m j n ~ 0.128 independent of T>\ this is expected from (11) recall- 
ing that x was found to be independent of D. The variance o 2 of the fitted 
distributions is linear in the noise strength V (see Figure 5(a)). This is also 
expected from (11) since x is independent of and a 2 is linear in V. 



The statistics for the time T between successive bursts are shown in Figure 4(c) 
for V = 1 x 1CT 7 . The pdf for T is fit very well by the Gaussian distribution 



p(T) 



exp 



-(T-T) 2 /(24) 



(14) 
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The variance a\ of the fitted distributions is linear in the noise strength T> (see 
Figure 5(b)), and T w 29.0 is independent of T>. It is interesting to compare 
these results to the results of [4,5] which considered the effect of noise on 
attracting, structurally stable heteroclinic cycles connecting finite amplitude 
fixed points. It was found that the distribution of times between successive 
visits near the fixed points on the cycle was not Gaussian; instead, it has 
a long, exponential tail for large times. There is a fundamental reason for 
the difference between those results and the results presented here: because 
the heteroclinic cycles considered in [4,5] connect finite amplitude solutions, 
the trajectory spends most its time in a small neighborhood these solutions. 
Indeed, this is what allowed a detailed analysis for the times between successive 
visits near the fixed points to be done. On the other hand, the heteroclinic 
cycles considered in this paper involve infinite amplitude solutions and are 
traced out in finite time [8]. More specifically, in the original time t defined 
by dt/dr = p, equations (8) become 

dx \ u x dp 
dt p ' dt 

Assuming the initial condition (xq, p) G S p at t = 0, these equations have the 
solution x(t) = x (p*/(p* - \ s t)) Xu/x % p(t) = -X s t + p*. The time to go from 
S p to Y, x is thus t = (p*/X s )(l - (x /x*) x °/ Xu ). For the limit x -> (i.e., 
as the heteroclinic cycle is approached), t — > p*/\ s , a finite time. A similar 
argument shows that the trajectory also returns from infinity in finite time 
(cf. [8]). This is why the time series in Figure 1(a) shows very rapid growth 
to and decay from large amplitude. We conclude that the time spent near 
solutions at infinity is not a dominant part of the total time between bursts. 
Thus, the results given in [4,5] are not relevant here. 

It is not possible to do detailed analytical work for the statistics of the times 
between bursts because this would require a more detailed understanding of 
how the trajectory behaves away from the solutions. However, the fol- 
lowing model for the flow between successive bursts gives some qualitative 
understanding. Consider the formal equation 

f t =v + v(t), (15) 

where v is constant, and n is a Gaussian white noise random process with 
(r}(t)) = and (r)(t)r}(t')) = 2VS(t-t'). If f = £ = at t = 0, the conditional 
probability distribution function P(£, t|£ ) for solutions of equation (15) obeys 
the Fokker-Planck equation 

d t P = -d^vP) + £%P, P = 5(0 for t = 0. 
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Fig. 6. (a,b) Trajectory for A = 0.1 and V = 0.01. (b) shows the detail near the 
solution at (p,9,(j)) = (0,7r/2, 0), and the arrows show the sense in which the 
trajectory is evolving. 



This has solution 



P(Z, *l& = 0) = : j=L= exp (-(£ - vtf/{AVt)) , 



i.e., a Gaussian centered at £ = vt and with variance 2Vt. The average time 
to reach a fixed value £ = £* is then T = £*/v, and the pdf for £ at this time is 
= P(£,T|£ = 0). We can re-express this in terms of a pdf for the time T 
at which £ = £* by using the deterministic approximation T^f-((-(*)/?). 
Then, using = |p(^)dT|, we obtain 

p(T) = exp {-v\T - Tf/iAVf)) . 

We make a connection to the distribution of times between successive bursts by 
thinking of £ as an arclength coordinate along a trajectory for equations (1,2). 
Because the trajectory does not spend a dominant part of its time in a neigh- 
borhood of the -Uqo and gpoo solutions, equation (15) is a reasonable model for 
the flow. Assume that the arclength traced out between successive bursts is 
roughly constant. The qualitative model (15) thus predicts a Gaussian distri- 
bution for the time between successive bursts, with the average time between 
bursts independent of and the variance linear in V. This is precisely what was 
found numerically. 

Figure 6 shows the projection of a numerical solution to equations (1,2) with 
the same parameter values as for the trajectories shown in Figure 1, except 
with D = 0.01. Unlike the result in Figure 1(b) for much smaller noise strength, 
here the underlying stable quasiperiodic solution is not even recognizable. The 
arrows in Figure 6(b) show that the trajectory can make a visit near a 
solution and depart either to the "left" or to the "right" . In the (p, x, y) 
coordinates introduced above (cf. equation (8)), the trajectory can intersect a 
surface of constant p for either positive or negative x values. Thus, sufficiently 
large noise can be responsible for the trajectory crossing the stable manifold 
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Fig. 7. Bursts for A = 0.1253 and (a) V = 0, (b) V = 1 x 1(T 7 . In (b), it is seen 
that the noise kicks the system out of the basin of attraction of the state shown in 
(a) and into the basin of attraction of a solution with constant p = 1/r. 

of the Moo solutions. This leads to the possibility that successive visits near 
infinity can either be to different or to the same Uoo solutions. In the convection 
system, successive visits near the same Uoo solution correspond to successive 
bursts occurring on the same side of the container [7]. Therefore, sufficiently 
large noise can affect the physical manifestation of the bursts by destroying 
purely blinking states. 

Reference [7] also identified a different type of stable quasiperiodic solution 
which is present in the absence of noise for A = 0.1253 (in the notation of [8], 
this is a u/v 1 solution). The trajectory for this solution makes successive visits 
always near the same solution (see Figure 7(a)). Physically, this corre- 
sponds to a winking state for the convection system because successive bursts 
occur at the same side of the container [7]. This solution coexists with stable, 
finite amplitude periodic solutions at (p,6,<f)) = (6.298, 1.753, 2.010 ± run), 
where m is an integer. Reference [7] argues that these stable, finite ampli- 
tude periodic solutions are more likely to be observed in experiments in which 
the Rayleigh number is ramped upwards, which helps to explain why winking 
states have apparently never been observed in binary fluid convection experi- 
ments. It is seen that even a very small amount of noise (T> — 1 x 10~ 7 ) can 
kick the system out of the basin of attraction of the winking state and into 
the basin of attraction of a stable periodic solution (see Figure 7(b)). Thus, 
even if a winking state could be established in an experiment, it would be 
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particularly sensitive to noise. 



4 Noise-induced bursting 

It will now be shown that related bursts may occur for equations (1,2) even 
in the absence of forced symmetry-breaking, that is, when noise is added to 
the normal form equations for the Hopf bifurcation with exact D4 symmetry. 
Specifically, consider equations (1,2) with the same parameters as considered 
in the previous section (A = 1 — 1.5i,B = —2.8 + 5i, C = 1 + i,u = 1), 
except with AA = Au) = 0. For these coefficient values in the absence of 
noise, a branch of unstable periodic solutions bifurcates subcritically from the 
trivial state, and two branches of unstable periodic orbits and one branch 
of globally attracting quasiperiodic solutions all bifurcate supercritically (see 
[15], also [7,8]). Infinite amplitude periodic and quasiperiodic solutions exist as 
counterparts of the finite amplitude solutions which bifurcate from the trivial 
state. In particular, and qp^ solutions are the infinite amplitude coun- 
terparts of the periodic solutions on the subcritical branch and quasiperiodic 
solutions on the supercritical branch, respectively [8]. If for some reason the 
trajectory comes close to the stable manifold of a Uoo solution, it will make 
an excursion to a neighborhood of that solution, then get kicked toward a 
qPoo solution which returns it to smaller amplitude. However, since the finite 
amplitude quasiperiodic solutions are globally attracting, heteroclinic cycles 
involving infinite amplitude solutions cannot form, so such a burst is a tran- 
sient phenomenon. A stable quasiperiodic solution is shown for A = 0.1, V = 
in Figure 8(a). For small noise the trajectory remains near the quasiperiodic 
solution (see Figure 8(b)). However, sufficiently large noise can repeatedly kick 
the system close to the stable manifold of a solution, leading to bursts of 
very large dynamic range (see Figure 8(c)). This is highly reminiscent of re- 
sults for the effect of resonant temporal forcing on the Hopf bifurcation with 
D4 symmetry [13]; there it was found that as the forcing amplitude increased, 
the attractor came closer to the stable manifolds of solutions, leading to 
bursting behavior. 

As a further demonstration of the effect of noise on a Hopf bifurcation with 
exact D4 symmetry, we close by considering an attracting infinite amplitude 
heteroclinic cycle. Specifically, consider the parameters A = 0.1 — 1.5i,B = 
-1 - 4i, C = 1 - 2i, A = 0.1, u = 1, AA = Auj = 0. Using the results of [15], 
it can be shown that there are heteroclinic connections between solutions 
within the £ subspace. Linearizing equations (4-6) about the solutions 
shows that the eigenvalue corresponding to perturbations in the p direction 
is —0.2, while the eigenvalues corresponding to perturbations within the £ 
subspace are —4.828 and 0.828. Since the first eigenvalue is negative and the 
sum of the latter two eigenvalues is also negative, the heteroclinic cycle con- 
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Fig. 8. Solutions for A = 0.1 and (a) V = 0, (b) V = 1 x 1(T 6 , (c) V = 1 x 1(T 4 
with AA = Aw = 0. Note the different scales for the time series plots. 



necting is attracting (cf. [15]). Figure 9 shows a numerically calculated 
trajectory for these parameter values and T> = 1 x 10 -5 ; the time step of 
integration is St = 1 x 1CT 6 . Much like the situations described in [3-5], noise 
causes the attracting heteroclinic cycle to become a "statistical limit cycle" 
in which switching between -u^ solutions occurs randomly in time but with a 
well-defined mean period. 
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0.03 U 




Fig. 9. The left figure shows the attracting infinite amplitude heteroclinic cycle 
connecting Uqo solutions at p = 0, and the "statistical limit cycle" solution for 
T> = 1 x 1CT 5 . The right figure shows that successive visits near Uqo solutions occurs 
randomly, but with a well-defined mean period. 

5 Conclusion 

The effect of additive white noise on a model for large aspect-ratio binary 
fluid convection has been considered. Particular attention has been paid to 
the effect of noise on bursting behavior present in the model. It was shown 
that even a very small amount of noise can have a very large effect on the 
amplitudes of successive bursts. This is because there will be a burst with 
larger amplitude if the noise kicks the system closer to the stable manifold of 
an infinite amplitude state, and a burst with smaller amplitude if the noise 
kicks it away. Numerical results and analytical arguments were given to un- 
derstand the statistical properties of bursts in the presence of noise. It was 
also demonstrated that large enough noise can affect the physical manifesta- 
tions of the bursts by destroying purely blinking and purely winking states. In 
fact, winking states are particularly sensitive to noise, which helps to explain 
why they have apparently never been observed in experiments. Finally, it was 
shown that related bursts can occur when noise is added to the normal form 
equations for the Hopf bifurcation with exact square symmetry. 

Reasoning similar to that given in [8] suggests that this type of bursting be- 
havior can persist when higher terms in equations (1,2) are retained, even in 
the presence of noise. Bursts will then be associated with visits near large but 
finite amplitude solutions. This is important because it implies that bursts 
similar to those described in this paper may be observed in real physical sys- 
tems undergoing a Hopf bifurcation with exact or weakly broken D 4 symmetry; 
the (unphysical) infinite amplitude solutions are not necessary. There is also a 
possibility that in appropriate cases the noise could be more important than 
higher order terms so their exact nature would be irrelevant. 

Other physical systems for which the results presented in the paper might be 
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relevant because their evolution equations are related to the Hopf bifurcation 
with D 4 symmetry include any system in a square domain undergoing an 
oscillatory instability [17] or displaying oscillatory patterns which are periodic 
on a square lattice [18], electrohydrodynamic convection in liquid crystals [19], 
lasers [20], spring-supported fluid-conveying tubes [21], the Faraday system in 
a square or nearly-square container [22], and dynamo theories of magnetic 
field generation in the Sun [10,23]. 
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